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■ We consider accretion of matter onto a low mass black hole surrounded by ionized medium. We 
show that, because of the higher mobility of protons than electrons, the black hole would acquire 
positive electric charge. If the black hole's mass is about or below 10 20 g, the electric field at the 

£\j ■ horizon can reach the critical value which leads to vacuum instability and electron-positron pair 

_ ' production by the Schwinger mechanism. Since the positrons are ejected by the emergent electric 

field, while electrons are back-captured, the black hole operates as an antimatter factory which 
' effectively converts protons into positrons. 

C<") ' PACS numbers: 98.62.Mw, 97.60.Lf 

■ Introduction - It is well-known fact that the difference between the mass of proton and electron can lead to 
predominant capture of protons by celestial bodies, which can make them electrically charged This process can be 
exemplified by the case of a cold object, such as a planet, surrounded by an atmosphere of positive ions and electrons 
in thermal equilibrium. Here, the escape velocity can be reached easier by electrons, which thus would leave the 
atmosphere. This process leads to the appearance of growing electrostatic field around the planet, making electron 

i escape more and more difficult. The equilibrium is eventually reached when the sums of all forces acting on ions 

and electrons are equal, which results in non-zero charge acquired by the planet. A similar mechanism can apply 
CNJ 1 to a radiating object. While protons and electrons are gravitationally attracted towards the radiating body with 
the same acceleration, the outgoing radiation acts differently on the two particle species. Indeed, as the interaction 
between the outgoing electromagnetic radiation and charged particles is described by the Thompson cross-section, 
the ej interaction dominates the P7 one by a huge factor of (m p /m e ) 2 ~ 3 • 10 6 , making proton-photon interaction 
£^ , effectively negligible. The stationary solution is again the one in which the astrophysical body is electrically charged, 
so the emergent electrostatic field makes protons and electrons fall with the same acceleration. In the standard case of 
radiative atmospheres the effect of the accumulated charge is negligible, so the electric field is not taken into account. 

In ref. @ , the authors discussed the case of a black hole (BH) which accreted spherically with luminosity close to 
the Eddington limit. If the BH mass is smaller than about 10 20 g, the electrostatic field of the stationary configuration 
at horizon exceeds the critical limit to produce electron-positron pairs. The crucial ingredient of the mechanism is 
that BH accretion has to be close to the Eddington limit, because the outgoing radiation has to inhibit the electron 
accretion. This is only possible if the density of the surrounding matter is huge, roughly above 10 protons/cm 3 rs 
1 g/cm 3 . In this paper we revisit this mechanism, taking the particle thermal distribution into account. We find that 
the process of proton-to-positron transformation can work even if the particle number density of the medium around 
the BH is much lower, close to realistic values. Interesting phenomenological implications are therefore possible. 

Schwinger mechanism at Schwarzschild horizon - Let us consider a BH of mass M surrounded by plasma 
of protons and electrons. In the simplest case of perfect spherical symmetry, the radial part of the equations of motion 
for the proton and electron fluids are respectively 

GnM aQ La lp a lp n^u:^ n v a pe P 

V P = 2^ + 12^ + ^2^ V P~— («*-*«)> (1) 
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GnM aQ La ie a~ e n~<J~ n e a pe P 
v e = 2 2 + ~a — T — —v e + ^ — (v p -v e ). (2) 

Here v p and v e are the proton and electron fluid velocities (that is the average velocities of p and e in the plasma 
which do not include the large chaotic thermal velocities), Q is the electric charge of the BH in proton charge units, 
a = e 2 /47r = 1/137, <Xy is the cross section of scattering of i on j, L is the luminosity in the comoving frame of the 
accretion flow, P is the momentum transfer in ep-scattering, n p and n e are the number densities of p and e around 
the BH, n 7 is the photon number density in the photon bath surrounding the BH and w 7 is the photon energy, which 
is roughly the momentum transfer in jyy— and e7-scattering. We neglect the angular momentum term because, as 
explained below, we are interested in the particles with low angular momentum. Eqs. |T]) and ^ indeed recover the 
Eddington limit for Q = and stationary flux ij — if the fourth term on the right hand side of the two equations, 
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that is, the term due to scattering on the thermal bath of photons, is neglected. Taking Q — 0, n p = n e and Vj = 0, 
we can solve the two equations for L, finding the usual Eddington luminosity, Le — 4:nG]\[Mm p /a e -y. 
In what follows, it turns out to be useful to define the dimensionless quantities 

K = — £ » 1 , q= — 2L and I = - ^— «1, (3) 

m e r g m p 4Trr g m p 

where r g = 2G^M is the gravitational radius of the BH. Since the cross sections of erf— and ^-scattering is inversely 
proportional respectively to m 2 and m 2 , the impact of the radiation pressure on acceleration for protons is suppressed 
with respect to electrons by K 3 . 

Let us consider particles inside the mean free path from the BH, r < Xj (j = p, e). If their velocity is small enough, 
i.e. Vj < Vj(r) = yjFj/r, where F p = r g (l — 2q — 21) and F e = r g (l + 2Kq — 21K 3 ), such particles would be captured 
by the BH. Even if we added the angular momentum term into eqs. fl} and (|2j). which is about J 2 /r 3 and has positive 
sign preventing accretion, the rather small velocities lead to a centrifugal to gravitational force ratio 

72 /„3 rv i p 

G^ = ft^7 < G^7 = 2 < 1 - 2 ''- 2, » < 4 » 

for protons and 2(1 + 2Kq — 21K 3 ) for electrons. Since this ratio for protons is always smaller than 2, the account 
of their angular momentum cannot significantly change the protons' accretion, while the electrons' accretion may 
noticeably slow down. This would increase the charging rate. Then, most of the captured particles are quickly 
swallowed by the BH. Here indeed the picture is quite different from the widely studied astrophysical accretion for 
large objects: charged particles can easily lose angular momentum by emitting classical electromagnetic radiation 
due to accelerating motion around the BH and probably even by Coulomb scattering, because of the higher particle 
number density in the vicinity of the BH. 

The particles propagate freely inside distance r, so we can assume thermal distribution and consider their motion 
under the influence of gravitational and Coulomb forces only. The number density of particles with velocity Vj < Vj 
is 

n c j (r) = C j (r) dvv 2 exp (-^— ). (5) 
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If the plasma were in equilibrium, Cj-(r) would be given by Cj(r) — Cj exp (G^Mraj /fT), where Cj is the value of 
Cj(r) far from the BH. However, in what follows we neglect the dependence of Cj(r) on r and we take Cj(r) = Cf 
the BH accretion is actually determined by the particle distribution far from the BH, where the chemical potential 
is small, while in the neigborhood of the BH, protons and electrons are not in equilibrium, since there is an intense 
ingoing flux to the BH. The quantity Cj can be determined by the normalization condition 

/W2T x3/2 
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Cj j ^exp(-— )=C j * r _) , (6) 



where rij is the particle number density at infinity and we can reasonably assume that n p — n e . 

The total number of particles with Vj < Vj inside the shell from r to r + dr is dnj(r) = ^Kr 2 n C j{r)dr. By definition 
of rij, drij(r) provides the number of particles at distance r which would be captured by the BH if they did not collide 
with other particles in the surrounding plasma. The account of collisions on particle propagation is described by the 
equation dn,j/dx = — Aj7i 7 , where Xj is the particle's mean free path in the medium. In the case of collisions with 
photons in the thermal bath with temperature T, the mean free path is given by 



A, = - (7) 

where <7j 7 = 87ra 2 /m|, n 7 = 0.24T 3 is the number density of photons and the last factor accounts for the necessity to 
make yj rrij / T collisions to transfer momentum of the order of the thermal momentum ~ yJrrijT . The latter factor is 
absent for low momentum particles. Anyhow its presence or absence does not significantly change the results. So, we 
multiply this quantity by exp(— r/Xj) to take into account the effects of the scattering in the plasma. To estimate the 
capture rate we need to divide this quantity by the time of the propagation from r to zero. This time can be found 
by the integration of the equation of motion fj = —Fj/(2r ). In principle, the result depends upon the initial velocity 
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and direction. However, an order of magnitude estimate for the time is t(r) — 2/3 r 3 / 2 Fj . The total number of 
captured particles (p or e) per unit time is 
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The main contribution to this integral comes from small x, where the exponent in the second term is equal to 1, so 

Nj = S^n.Ff ( £jL) In ^ . (9) 
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The temperature of protons may be smaller than the temperature of electrons. The latter is usually equal to the 
temperature of the photon bath. The mean free paths Xj of e and p are also different, but their impact on the result 
is weak, because Nj depends on Xj logarithmically. 

The equilibrium state would be reached if N p = N e . Keeping in mind that F p « 1 + 2q and F e ss 1 + 2Kq, we find 
that the equilibrium takes place for 

(1 + 2qf{m p /2T p fl 2 HX p /r g ) « (1 + 2Kq) 2 (m e /2T e fl 2 ln(A e /r 9 ) . (10) 

Let us assume Kq ^ 1 ^S> q. In this case we get 

4K 2 q 2 * K^ 2 f^f 2 . (11) 

\T P J ln{Xe/r g ) 

Since X p /X e ~ K 5 / 2 (see eq. (JT))), at the temperature T = 1 keV one finds 

ln(A p /r g ) _ 1 | ln(A p /A e ) _ 1 ? ^ 
ln(A e /r 9 ) ln(A e /r 9 ) 

Thus Kq ~ 100 and g ~ 0.1 and the assumption Kq ^> 1 » q is consistent. The accumulated electric charge at 
equilibrium is 



~ (Te/Tp) v Rv^) (13) 

and, for T p = T e , we obtain aQ « 0.10 m p r g . The electric charge Q, and hence the efficiency of our mechanism, 
strongly depends upon the temperatures of protons and electrons. We expect that they are different for different 
astrophysical environments. The case T p = T e = T 7 is the most reasonable, because the coupling of electrons and 
protons to the thermal bath of the surrounding photons is strong enough to establish equality of their temperatures. 
If it were not so, one may argue that T p is larger than T e because the accreting protons have virial temperature, 
but electrons cool efficiently, so that T p /T e 1 and Q is much smaller than our estimate. That does not occur 
here, because the accretion onto small BHs is quite different from the one onto the Solar-mass BHs. The Eddington 
accretion for small BHs is typically several orders of magnitude larger than the Bondi accretion, because the former 
is proportional to the BH mass, while the latter - to the square of the BH mass. For example, in the case of 10 20 g 
BH inside a cloud with density 10 14 cm~ 3 , the Eddington to Bondi accretion rate is about 10~ 10 . We should thus 
expect that there are not many collisions between particles of the accreting matter and the gravitational energy is 
not released efficiently, but advected into the horizon where it is lost. The efficiency parameter 77 is surely very low 
and T e is close to T p . Moreover, it is even more probable that there is no accretion disk at all: the gravitational 
radius of a 10 20 g BH is about 10" 8 cm, i.e. the size of an atom. At such small scales, quantum effects are important 
and they should prevent the formation of a "regular" disk. On the other hand, as already said, most of the captured 
particles are quickly swallowed by the BH because they emit radiation to their orbit acceleration and thus lose angular 
momentum efficiently. 

If the electric field at the horizon is strong enough, there is the possibility of particle pair production by vacuum 
breakdown (or Klein instability) [3J . If the BH radius is much larger than the Compton wavelength of the electron 
l/m e — 4 • 10~ n cm, i.e. the BH mass is M 3> 2 • 10 17 g, we can assume that the electrostatic field near horizon does 
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not change with the distance and the electron-positron production proceeds according to the Schwinger result 
The pair production probability per unit time and volume is 



w = -4- — 
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where E c = rn^/e. Even if E ^ implies W ^ 0, the production of particles in a uniform electrostatic field is efficient 
only if E is close to the critical value E c . In our case, the equilibrium electric field at horizon is 
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and E c /E = 0.9 (M/10 20 g). The accretion rate of protons © is 



p ~ U 10 10 /cm 3 
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and since the value of the electric field at horizon is close to E c , the e + e~ pair production is fast and the rate is at 
the level of the proton accretion rate, that is eq. (Til)]) . Then, the created positron would run away from the BH due to 
Coulomb repulsion, while electrons are predominantly back-captured: the BH works actually as an efficient antimatter 
factory, converting protons into positrons. The energy of positrons at infinity is about aQ/r g 0.15 m p ~ 140 MeV. 

DISCUSSION - Small BHs with the mass at the level of 10 20 g surely cannot be produced by the stellar collapse, 
but might be created in the early Universe, through many possible mechanisms [H, [6j]. Then, if their mass is larger 
than 5 • 10 14 g, their lifetime exceeds the age of the Universe and they could still live somewhere today. The proposed 
mechanism of positron production works if the primordial BH is surrounded by a ionized medium with the density of 
roughly 10 10 particles/cm 3 , which is about fourteen orders of magnitude smaller than the density required in ref. (2j. 
For example, the mechanism can operate in the neighborhood of the super-massive BH in the Galactic Center. Here 
one finds an ionized medium with temperature of roughly 1 keV and particle number density of the order of 10 s cm~ 3 
in a sphere of radius ~ 10 16 cm around the central BH [7|. A more favorable environment is the hot and dense 
atmosphere around accreting compact stars: in the accretion disk, the temperature of the plasma can reach 10 keV 
and the particle number density can be as high as 10 20 cm -3 8]. Primordial BHs may be abundant around "old 
compact objects", e.g. the final product of the first stars, because captured in the early Universe @. They can still 
live today in the accretion disk of these objects because they cannot be swallowed quickly by their hosts like ordinary 
matter: they do not lose energy easily. The main mechanism of energy loss for small BHs in a plasma is via dynamical 
friction and the phenomenon is discussed in ref. [§]. One can see that, for BHs with a mass around 10 20 g and for 
reasonable values of density and temperature of the plasma, the time necessary to fall into the central object exceeds 
the Universe age. 

In principle there are three observable phenomena induced by the process of proton-to-positron transformation. 

i) Observation of MeV positrons in cosmic rays: such positrons would be hard to observe, because they would 
quickly stop and annihilate in the high density environment. Indeed the mean free path of 100 MeV positrons due to 
scattering on the photon bath with temperature T is 
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So they stop quite fast and annihilate with electrons. Also, the cross section cy should be smaller than the Thompson 
one, but for 100 MeV positrons and 1 keV photons the correction is of order one and does not change our conclusions. 

ii) Observation of the 0.511 MeV line from e + e~ annihilation. Here the problem is that the dense medium, 
where the small BHs should be to efficiently transform protons into positrons, would disperse this line to a continuos 
background. Still the line might survive if the effective size of the dense region around the BHs is, at least in one 
direction, smaller than the photon mean free path. The latter can be estimated as 
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for 0.5 MeV photons due to scattering on the background electrons. So, the possibility of observing the line depends 
crucially on n e . It is natural to wonder whether the mechanism can explain the observed 511 keV line from the 
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Galactic Center [1(3], where the expected annihilation rate is at the level of 10 43 s 1 and several different mechanisms 
have been already proposed [ll[ . From eq. (fT6|) , we see that the number of BHs should be 

Since the total mass of non-baryonic matter in the Bulge cannot be larger than 10 10 Mq 12], the small BHs have to 
live in regions with n p > 10 14 cm -3 . It is difficult to make concrete predictions, as the number density of BHs and the 
temperature and particle number density of the environment are unknown. The upper bound on Nbh is given by the 
condition that their total mass cannot be larger than the maximum mass of dark matter (there are indeed no other 
constraints for BHs with masses in the range 10 17 — 10 27 g @). Their distribution would be surely inhomogeneous 
and it is not unreasonable that they can be accumulated around compact bodies in the center of the Galaxy. 

Hi) A more realistic possibility is likely the observation of the total gamma ray luminosity created by the annihilated 
100 MeV positrons. According to eq. (fT75|) . each BH emits about 10 16 erg/s for n p = 10 14 cm~ 3 . If a total amount of 
10 — 100 Mq small BHs are accumulated into a high density cloud/atmosphere, they may be observable as a luminous 
object with the power of about 10 30 — 10 31 erg/s. 

Conclusion - If primordial light black holes surrounded by a ionized medium exist in the Galaxy, they could 
operate as an efficient antimatter factories converting accreting protons into positrons. The charging of BH by the 
predominant influx of protons in comparison to electrons is achieved because of the much larger proton mass and 
correspondingly much larger proton mobility in the surrounding medium. If the black hole is small, the electrostatic 
field at the horizon can exceed the critical value of the vacuum stability and electron-positron pair production by 
Schwinger mechanism becomes efficient. Primordial black holes with the mass in the range 10 15 — 10 20 g might live 
today in the Universe, injecting positrons. If these black holes live in an ionized high density medium, such as around 
the super-massive black hole in the Galactic Center or in the atmosphere of violent stars, the proton conversion 
into energetic positrons would enrich cosmic rays by positrons with energeis in the interval 1-100 MeV and gamma 
radiation of similar in energies, which possibly makes the effect experimentally observable. 
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